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Recent developments in the evaluation of two-loop pseudo-observables and observables are briefly reviewed. 



1. Introduction 

The complete strategy to derive theoretical pre- 
dictions for pseudo-observables (PO) and observ- 
ables (0) up to two-loop accuracy is based on the 
following steps: generation and manipulation of 
diagrams, rcnormalization, semi- numerical eval- 
uation of diagrams. We perform the first step 
with the help of the code QraphShot (FORM 
3.1) [ P while the last is performed using the 
code CoopBack (FORTRAN 95) which makes 
extensive use of array handling, assignment over- 
loading, vector/recursive functions. 

2. Renormalization and Unitarity 

A keyword in renormalization is counterterm; 
they are not strictly needed but are, neverthe- 
less, very useful when dealing with overlapping 
divergencies. QraphShot generates all countert- 
crms needed in the standard model and produces 
ultraviolet (UV) finite Green functions. Another 
important keyword is skeleton expansion, mean- 
ing that the relevant objects in perturbation the- 
ory are dressed propagators. With their help we 
want to construct a (finite) rcnormalization pro- 
cedure where a) the renormalized parameters are 
real, b) finite renormalization is the result of a 
consistently truncated solution of renormalization 
equations (RE), c) complex poles arise naturally 
after dressing the propagators, but cutting equa- 
tions remain valid to all orders. For an unstable 
particle V define 



(1) 



where S is the V (skeleton) self-energy. Cutting- 
equations and unitarity of the S -matrix can be 
proven: one uses two-loop A in tree diagrams, 
one-loop A in one-loop diagrams and tree prop- 
agators in two-loop diagrams. The proof is due 
to Veltman [|3]: the crucial observation is that A 
satisfies the Kallen - Lehmann representation, 



At(p 2 ) = 6(po) A v (p 2 ) 2Im£ vv (p 2 ) 



(2) 



while, for a stable particle s, the pole term shows 
up as 



A+(p 2 ) =% ) [A s (p 2 )\ 2Im£ ss (p 
+ 2inS(p 2 + m 2 ). 



(3) 



One then expresses ImE vv m terms of cut self- 
energies, repeats ad libidum and derives that con- 
tributions from cut lines come from stable parti- 
cles only. Consider a toy model with 



L Snt = |$(x) <t> 2 {x), 



(4) 



1 A v S vv 



and where 4> is unstable. We define A$ and A^ 
according to Eq.QJ. An example of the skele- 
ton expansion for self-energies is given in Fig. 1: 
Im 7^ only due to the 3— particle cut of di- 
agram b) of Fig. 1 and only diagrams a) and c) 
are retained in the expansion; in a) we use A$, 
at one-loop accuracy. 

In a gauge theory, however, there is a clash be- 
tween resummation and gauge invariance; usually 
only the complex pole is rcsummed [ 0] . There- 
fore a one-loop self-energy with A$ at one-loop 
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accuracy is equivalent to the 3 diagrams of Fig. 2 
computed with A$(s M ), where 

g 2 

z p = T^ B o (~ s Af ; m , m) . (5) 

lb 7T Z 

The interplay of gauge parameter independence, 
Ward - Slavnov - Taylor identities and unitar- 
ity (which is naturally satisfied in the framework 
of dressed propagators) requires a more detailed 
analysis, beyond the scope of this paper. 

To solve REs we need an input parameter set 
(IPS) including some notion of M exp . In the past 
on-sell PO have been derived by fitting lineshapes 
from experiments [Ej but we can use an on-shell 
M os only at one-loop. Beyond this order the 
correct treatment requires introducing complex 
poles [ |B] , s v = fi% — i "f v fi v . 

If we want to use available data a transforma- 
tion is therefore needed: define pole PO through 
-0 = arctanr os /Af os : M P = M os cosi/j and 
T P = r os sin^. Once again, a change of strategy 
is needed since RE change their structure at two 
loops. It is a new perspective: at one loop one 
considers M os as an input parameter indepen- 
dent of s P and derive s P . At two loop REs are 
written for real p R and solved in terms of (among 
other things) experimental s P . 

In the framework of an order-by-order renor- 
malization, M R is a real solutions of truncated 
REs consistently with cutting-equations and uni- 
tarity. 

2.1. Complex poles 

Here we give an example of the old fashioned 
one-loop technique: on-shell masses are input 
and, having the Higgs boson in mind, we derive 
s h — — t/i|[7 H ; numerical results for s H are 
shown in Tab. 2] 

The current fashion instead is to extract s^f p = 
[i 2 H — i [i H j H from data and to derive s^ 1 = M§ — 
iM H T H . To see how it works we point out the 
difference with previous two-loop calculations [d 
where one starts from 

s v = to 2 - T, vv (s v , m 2 , . . . ) , (6) 

and derives the complex pole in terms of the 
(bare) renormalized mass 

s v = m 2 — E^l (to 2 , m 2 ,...)+.. . (7) 



Table 1 

One-loop s H (M° s ) 



M° s [GeV] 


120 


300 


Mh [GeV] 


119.96 


299.74 


7h [GeV] 


7.00 x 10~ 3 


7.90 



To improve the quality of the result we consider 
the relation 

s v = to 2 - Y, vv (s v , to 2 , {p} , . . . ) (8) 

where with {p} we denote additional, renormal- 
ized, parameters. We then solve REs and obtain 
to 2 and {p}, i.e. 

m 2 , M = Re/(s Vl , s V2 , . . .) , (9) 

where s Vi are experimental complex poles. Note 
that we never expand functions depending on 
complex POs which means that we actually com- 
pute two-loop diagrams on the second Riemann 
sheet. Furthermore, in our scheme the solution 
of REs is used at the Born level in two-loop di- 
agrams, one-loop in one-loop diagrams, two-loop 
in tree diagrams. If V £ {Vi , V% , ■ ■ ■ } we have 
a genuine prediction, otherwise we have a con- 
sistency relation for loop corrections. Numerical 
results for the two-loop s^ 1 (s^, xp ) are shown in 
Tab.0 

Table 2 

Two-loop sf(s%* p ). All entries in GeV. 





1h 


M H 


r H 


300 


4 


299.96 


8.374 


300 


12 


299.87 


8.376 


500 


40 


500.17 


63.37 


500 


80 


500.42 


63.34 
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Table 3 

Percentage two-loop corrections in RE Ea. (|10f> . 



M° s [GcV] 


150 


300 


500 


G F rf v £<f> 

2tt 2 TUT 

°G 


18.29% 


8.89% 


-24.62% 



2.2. Numbers & renormalization 

Consider one of the REs, e.g. 

2 

(1 + A 5 ), A 5 = <5 G + A.g s 



G 



g 



(10) 



V2 8M 2 

relating g to the Fermi constant G F ; Ag s is the 
W self-energy part. A solution in perturbation 
theory starts with 



g 2 = 8G F vl 

1 
2 



1 + CW % 



Resl v | v (s w ) — ^^^(O) 



(11) 



According to an old result 6 G is UV/IR finite; 

(2) 

we can add that S G is finite after one-loop renor- 
malization. Furthermore, we define a process in- 
dependent Fermi coupling at the two-loop level 

G = G F {l-5P^4 L +\2(5o 1) f 
I 2 -k z L 

G F fly 



— s {1} <?« 

,.2 ° G °9 

H'w 



•,ir 2 J J 



(12) 



Therefore, starting from a RE like 
X = x (1 + ai a; + a 2 x 2 ), 



X 



_ G F fly 

2 IT 2 
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6^ + Sm, 



5f 



a 2 = ^ > 
we obtain the following solution: 

x = X + X 2 (b 1 +b 2 X), S {n) = 



(13) 



(14) 



The LO/NLO/NNLO terms are X, biX 2 and 
b 2 X 3 ; results are shown in Tabs. where one 
can see that perturbation theory becomes ques- 
tionable beyond 350 GeV. 



Table 4 

Percentage two- loop corrections in RE Eg. (11411 . 



M° s [GeV] 


150 


250 


350 


NLO/LO (%) 


+3.31 


-2.30 


-7.85 


h 


+ 12.28 


-8.51 


-29.07 


b 2 X 


+0.25 


-1.38 


-9.26 


NNLO/NLO (%) 


+2.06 


+ 16.16 


+31.85 



3. Running of a 

The role played by the running of a has been 
crucial in the development of precision tests of the 
SM. Underneath this concept there is the popular 
wisdom that universal corrections are the impor- 
tant ingredient while non-universal ones should 
be made as small as possible; therefore, univer- 
sal corrections should be linked to a set of POs 
and data should be presented in the language of 
POs which, in turn, is connected with resum- 
mation, against gauge invariance. Admittedly, 
around M z it has been easy to perform a dis- 
crimination, relevant vs. irrelevant terms, paying 
a little price to gauge invariance. Well above this 
scale the situation is drastically different. Thus, 
the natural question is about the definition of the 
running of a at an arbitrary scale. One (fuzzy) 
idea is to import from QCD the concept of MS 
couplings and to express theoretical predictions 
through MS couplings [ |H] ■ This idea is open 
for criticism: although the MS parameter seems 
unambiguous it violates decoupling. 

There is another, well-known, solution: do the 
calculation in the gauge, select a £ - inde- 
pendent part of self-energies, perform resumma- 
tion while leaving the rest to ensure independence 
when combined with vertices and boxes. The ob- 
vious criticism is: it violates uniqueness; however, 
it is only a matter of conventions. 

Ingredients for ajjg are: a bosonic part, a 
fcrmionic part with 3 lepton generations, a per- 
turbative quark contribution (top or diagrams 
where light quarks are coupled internally to mas- 
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sive vector bosons) and a non-perturbative one 
with diagrams where a light quark couples to a 
photon (related to Aa;j ad (Mj)) . We define 

1 

47T 

Alternatively, we consider £ = 1 and define 



MS" 



a T zJ a ) = a-L-—TL™(p) 



(15) 



a(s) 



1 + Aa(s) 



(16) 



Numerical results are shown in Tab. 

4. Virtual infrared corrections 

We have been able to prove [0] that two-loop 
vertices have an integral representation 



dC k ({x}) 



A 



In 1 



or Li„ 



(17) 



where v4,j3 are multivariate polynomials in {a;}, 
the Feynman parameters. Two - loop diagrams 
are always reducible to combinations of integrals 
of this type where the usual monomials that 
appear in the integral representation of Nielsen 
- Goncharov generalized polylogarithms are re- 
placed by multivariate polynomials of arbitrary 
degree. 

Fig. 3 Example of IR divergent two-loop vertex. 



V K 



1>2 



M 



-P 




Pi 

Our method is fully multi-scale, it allows for a 
classification of infrared divergent configurations, 
for the evaluation of IR residues and IR finite 
parts and is also suitable for collinear regions. To 
reach these objectives we had to extend Berstein - 
Sato - Tkachov functional relations [E] to higher 
order transcendental functions. Results for the 
diagram in Fig. 3 are shown in Tab. 
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Table 6 

Comparison with the results of Davydychev - 
Kalmykov [E3j- Only the infrared finite part is 
shown, in units of 10 -8 GeV~ 4 . 







RcV , K 


Im V ;k 


Our 
DK 


400 


5.1343(1) 
5.13445 


1.94009(8) 
1.94008 


Our 
DK 


300 


5.68801 
5.68801 


-1.61218 
-1.61218 


Our 
DK 


200 


9.36340 
9.36340 


-2.84232 
-2.84232 


Our 
DK 


100 


29.4726 
29.4726 


-9.74218 
-9.74218 
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Figure 1. Example of skeleton expansion. 



a) skeleton 



b) S insertion 



<X> 



c) skeleton 

Figure 2. Rearranging perturbation theory in the 
presence of complex poles. 
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